A diakoptic approach to quantum computation 

Giuseppe Castagnoli*, Dalida Monti"*" 
In the diakoptic approach, mechanisms are divided into simpler parts 

interconnected in some standard way (say by a "mechanical connection"). We 

explore the possibility of applying this approach to quantum mechanisms: the 

specialties of the quantum domain seem to yield a richer result. First parts are 

made independent of each other by assuming that connections are removed. 

The overall state would thus become a superposition of tensor products of the 

Q ■ eigenstates of the independent parts. Connections are restored by projecting 



On 



(N 



off all the tensor products which violate them. This would be performed by 

particle statistics, under a special interpretation thereof. The NP-complete 
> 

problem of testing the satisfiability of a Boolean network is approached in 

p,i this way. The diakoptic approach appears to be able of taming the quantum 

t^^ ■ whole without clipping its richness. 

PACS: 89.70.+C, 89.80.+h. 






I. DEFINITION OF QUANTUM MECHANICAL CONNECTION 
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In (classical) applied mechanics, the diakoptic (dissectionistic) approach is exemplified 
by the notion of mechanical Connection. Connections divide the whole into simpler parts 
and reconstruct it — they introduce a "divide and conquer" strategy. In fig. 1(a), a crank- 
shaft is the Connection which imposes an invertible function between the positions of parts 
r and s (discretized as and 1, then the function is the Boolean NOT). 

Things can be more difficult in quantum mechanics, since the Connection Hamiltonian 
may not commute with the parts Hamiltonians. This difficulty is avoided by implement- 
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ing each Connection through a form of constructive and destructive interference, assumedly 
related to particle statistics. By applying reverse engineering, the Connection is first intro- 
duced as a mathematical feature that would be nice-to-have in quantum mechanisms. Then 
we ask ourselves whether that feature can be physical. 

Let us consider the mechanism of fig. 1 from a quantum mechanical perspective. The 
Connection should establish a constraint between two otherwise independent quantum parts 
r and s, with eigenstates |0) r , |l) r and |0) s , \l) s (fig. lb). The Connection state should have 
the form 

W) = a |0) r |l) s + P |l) r |0), , with \a\ 2 + \(3\ 2 = 1. 

We should note that the eigenvalues of each tensor product satisfy the Boolean NOT — 
the constraint; \<p) is free to "move" in a two-dimensional Hilbert space: this gives the one 
degree of freedom required from a Connection. 

Let us assume that the Connection is temporarily removed. The generic state of the in- 
dependent parts are |\&) r = a r |0) r + /3 r |l) r , |\P) S = a s \0) s + j3 s \l) s . The whole unentangled 
state in the Hilbert space of the two qubits H w is 

l*> = «o |0) r |0) s + ai |0) r |l) s + a 2 |l) r |0) s + a 3 |l) r |l) s , 

with a = a r a s , etc. The Connection is restored by projecting |\P) on the "symmetric" 
subspace H s = span {|0) r |l) s ,|l) r |0) s }. Let us define the projector (or "symmetry") A rs by: 

^ s |o) r |i) s =|o) r |i) s , A-.|i> r |o>, = |i> r |o>,, 
^ s |o) r |o) s = A-.|i> r |i>, = o. 

The A rs projection of |\1/) is the normalized vector of H s closest to it. This is obtained (in 
a peculiar way whose motivation will be clarified) by submitting a free normalized vector \<p) 
of H w (whose amplitudes on the basis vectors of H w are free and independent variables up 
to normalization) to the mathematically simultaneous conditions: (i) A rs \tp) = \tp) , and (ii) 
the distance between the vector before projection |\&) and that after projection \<p) should 



be minimum; in equivalent terms ||(\l/ \<p)\\ should be maximum. This yields the usual result 
\ip) = k («! |0) r |l) s + «2 |l) r |0) s ) , where k is a renormalization factor, namely an allowed 
Connection state. The Connection will perform by operating on the parts under continuous 
A rs projection on H s of the whole. This will turn out to be the nice-to-have mathematical 
feature. 

II. A DIAKOPTIC INTERPRETATION OF PARTICLE STATISTICS 

To give an introductory example, let us introduce a sort of Connection simply related to 
particle statistics. Let 1 and 2 be two free, identical and non-interacting spin 1/2 particles. 
At a given time, their overall spatial wave function is a symmetrical/antisymmetrical linear 
combination of the spatial wave functions of the two free particles, (x\ and x 2 are the particles 
spatial coordinates): 

the + (— ) sign goes with the spin singlet (triplet) state (normalization is disregarded). It 
can be seen that ||\P (xi,X2)|| — cos 2 kx for the singlet state, and ||\P {x\X 2 ) \\ = sin 2 kx for 
the triplet state, where x = x\ — x 2 , k = kA — hs- Thus close (separated) particles are 
more likely to be found in a singlet (triplet) state. There is a sort of Connection inducing 
a correlation between the mutual distance of the two particles and the character of their 
spin state: in principle, by operating on the distance, the character of the spin state is 
(probabilistically) changed. Noticeably, this kind of Connection would fall apart if the two 
particles were not identical. 

The Connection of Section I is a different case. It still relies on particle statistics, but 
under a special interpretation thereof. A particle statistics symmetry should be seen as 
the result of continuous projection of the system state on a given symmetrical subspace. 
This interpretation can be exemplified by considering a pair of identical bosons labeled 1 
and 2; 512 = \ (1 + Pvi) is the usual symmetrization projector; 0/1 stand for, say, horizon- 
tal/vertical polarization. The symmetry Syi |^) = |^) is satisfied in 



H t = span \0) 1 |0> 2 , |l) x |1) 2 , -= (|0> 1 |1) 2 + |1> 1 |0) 2 ) 
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There is a common didactic way of introducing this kind of symmetry. First, statistics 
is disregarded and the particles are assumed to be independent of each other. Let their 
unentangled state at some given time be |\l/ (£)) = «o |0) 1 |0) 2 + «i |0) 1 |1) 2 + «2 1 1) x |0) 2 + 
ct3 |1) 1 |1) 2 . Second, statistics is recovered by symmetrizing |\P (t)) , namely by projecting 
it on H t . We take this didactic procedure seriously: particle statistics is interpreted as the 
result of projection of the system state on a predetermined Hilbert subspace, the one which 
satisfies the symmetry. This amounts to considering the equation 

\/t:S 12 \*(t)) = \V(t)), (1) 

as a constraint applied to |\P (£)). When a particle statistics symmetry is an initial condition 
conserved as a constant of motion, this constraint is redundant. However, the notion of 
Connection will be related to particle statistics by means of a counterfactual reasoning 
based on eq. (1). The idea is that |\& (£)), symmetrical at time t, could be pushed out of 
symmetry at time t + dt; but in this case eq. (1) would "immediately" project it on H t back 
again. Particle statistics would operate like a watch-dog effect internal to the endosystem, or 
like destructive and constructive interference, by killing the amplitudes of those eigenstates 
of |\P (t + dt)) which violate the symmetry, and reinforcing the other amplitudes through re- 
normalization. This can also be seen as a continuous form of partial state vector reduction 
on a symmetrical subspace. 

To see why |\I/ (t)) could be "pushed out of symmetry", we must consider the system 
defined in Section I and the A rs symmetrization projector. In a first step, A rs projection is 
disregarded while parts r and s are assumed to be independent of each other. An operation 
on part r could well push the overall state |\P (£)) out of symmetry, but in a second step this 
is prevented by the continuous projection of \^f (t)) on H s : Vt : A rs |\P (£)) = |\P (t)) . 

We should note that this projection (or, if one prefers, partial state vector reduction on 
a predetermined subspace) will in general alter the entanglement between the parts r and 



s, thus the coherence elements of p r (t) (part r density matrix). However, it does not alter 
the diagonal of p r (£); this is determined by the operation performed on part r, namely it is 
a constraint that should be satisfied by projection. 

III. BEHAVIOUR OF THE QUANTUM MECHANICAL CONNECTION 

We go back to the Connection r,s defined in Section I, and consider an operation 
performed on just one qubit, say, r. Let this be the continuous rotation cosip\0) r (0| r — 
sin p |0) r (l| r + sirup |l) r (0| r + cosy? |l) r (l| r , with <p — cot and time t ranging from to ^-. 
We shall examine the effect of applying Q r (ip) to qubit r, 

Pr (t) = Q r {uot) pr (0) Ql (wt) , (2) 

under continuous A rs projection of the overall state. 

Let the Connection initial state be the "symmetrical" state (whose tensor products satisfy 
symmetry A rs ): 

|*(O))=cos0|O) r |l), + sin0|l) r |O),. (3) 

The successive states are obtained by submitting a free normalized vector |\P(£)) of the 
Hilbert space H w (Section I) to the following mathematically simultaneous conditions, 
for all t or p: 

i) A rs \y(t)) = \y(t)); 

ii) p r (t) = Tr s (|tf (£)> (V (0|) = cos 2 (0 + p) |0) r (0| r + sin 2 (7? + p) |l) r (l| r ; Tr s means 
partial trace over s. Under condition (i), |\P (t)) has always the form a (t) |0) r |l) s + 
j3 (t) \l) r |0) s , therefore p r (t) is a diagonal matrix: its coherent elements have been 
killed by A rs projection, or reduction; 

iii) the distance between the vectors before and after projection must be minimum. In 
the case of continuous projection, ||(\1/ (t)\ ^f (t + At))|| must be maximized orderly for 



t = 0,t = At,t = 2 At, ..., t = NAt, where At = $£; then the limit for JV -> oo 
should be taken (however, maximization ordering turns out to be irrelevant here). 

(i) and (ii) yield |\P (£)) = cos (■# + <p) |0) r |l) s + e l<5 sin (■# + </?) |l) r \0) s , with <5 unconstrained, 
as can be checked; condition (iii), given the initial state (3), sets 5 = 0, yielding to the unitary 
evolution: 



|* (£)> = cos (tf + <p) \0) r |l) s + sin (tf + ^) |l) r |0) ; 



(4) 



This makes a "good" Connection. The rotation of qubit r is identically transmitted to 
the other qubit s. In fact 



Tr r (|* (i)> (* (t)\) = p s (t) = sin 2 (0 + <p) |0) s (0| a + cos 2 (t? + ^) |1>, (1|, . 



(5) 



Of course eigenvalues and 1 are interchanged: one qubit is the NOT of the other. Notice- 
ably, by simultaneously rotating also the other extremity s of the Connection by the same 
amount, the same result (4) is obtained. This means adding eq. (5) as a condition, but 
this is redundant with respect to (i) and (ii), it was derived from (i) and (ii). Whereas, 
two different rotations of the two Connection extremities give an impossible mathematical 
system; this resembles a rigid classical Connection. 

It should be noted that a rotation ip of qubit (part) r under A rs projection, is equivalent 



to applying the unitary operator Q (ip) to |\& (£)): 

cosip sirup 



Q(v>) 



— sin ip cos (p 

cos ip — sin ip 



V 



,with|0) r |l>. 
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V 



sirup cos <p 

Q (ip) brings from |\& (0)) (3) to |\& (t)) (4) without ever violating A rs . We have thus ascer- 



v 1 / 



\ a i 



tained a peculiar fact. Our operation on a part, blind to its effect on the whole, performed 
together with continuous A rs projection, generates a unitary transformation which is, so to 
speak, wise to the whole state, to how it should be transformed without violating A rs . Of 
course A rs ends up commuting with the resulting overall unitary propagator (shaped by it). 
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IV. QUANTUM COMPUTATION NETWORKS 

Let us consider the reversible Boolean network of fig. 2(a), fully deployed in space — 
time is orthogonal to the network lay-out. This is different from sequential computation, 
where the Boolean network appears in the computation space-time diagram. 

Nodes t, u, v and r make the input and the output of a controlled NOT; r and s belong 
to a Connection. This c-NOT is made up of four coexisting qubits, and has four eigenstates 
which map the gate Boolean relation and constitutes the basis of 

H g = span{\0) t |0>„ |0)„ |0) r , |0) t |l) u |0)„ |l) r , |l) t |0) u |1>„ |l) r , |l) t |l) u |1>„ |0)J. 

Model Hamiltonians of such gates are given in [1,2]; this is different from time- sequential 
gates where the input and output are successive states of the same register J 3 ' 4,5 ' 6 ! 

The satisfiability problem is stated by constraining part of the input and part of the 
output, and asking whether this network admits a solution. Let u — 1 and s = 1 be such 
constraints, u — 1 (s — 1) propagates a conditional logical implication from left to right 
(right to left). Logical implication is conditioned by the values of the unconstrained part of 
the input (output). To have a solution, the two propagations must be matched, i.e. they 
must generate a univocal set of values on all the nodes of the network. Finding whether 
the network admits at least one match (one solution) is an NP-complete problem. Possible 
collisions (mismatch) between the two propagations will be both overcome and reconciled 
by the Connection. 

Let us assume that the network has just one solution (which is the case here: t — 1, 
u = 1, r = 0, v = 1, s = 1). The procedure to find it is as follows (this will hold for a 
generic network, thus we will think of many gates and Connections — in fig. 2b each wire 
is a Connection). The output constraint is removed while an arbitrary value, here t — 0, is 
assigned to the unconstrained part of the input. The logical propagation of this input toward 
the output yields t — 0, u = 1, r — 1, s = (v = t will be disregarded). This computation 
is performed off line in polynomial time. It serves to specify the initial state in which the 
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network must be prepared: |\& (0)) = |0) t |1) M |l) r |0) s . This state satisfies the gate/s and 
the Connection/s, but qubit s is in |0) (0| rather than |1) (1| (the output constraint). 
It will be continuously rotated from |0) (0| to |1) (1| under A rs projection, while keeping 
p u = \l) u (l\ u fixed. This transformation operates on the network Hilbert space H n ; here 
H n = H g ® H s where H s = span {|0) s , |l) s }. Note that all states of H n natively satisfy the 
gate/s, not necessarily the Connection/s. 

At any time t, the state of the network is obtained by submitting a free normalized state 
|\& (t)) of H n to the conditions: 
for all t: 

i) 4, |tt(t)) = |*(f)>; 

ii) Tr tiU , r (|*(t))(W(t)|) = p s {t) = cos 2 v?|0) s (0| s + sinV|l) s (l|, with <p = ut and t 
going from to j-; 

iii) Tr t , r ,s (1^ (0) (^ WD = P« (0) = l^)« (^Li ^ n a generic network there might be more 
conditions of this kind; 

iv) the distance between the vectors before and after projection should be minimum, as 
specified in Section III. 

This yields: |tf (£)) = cos <p \0) t \l) u |l) r |0) a +e lS smip\l) t \l) u |0) r |l) s , with <p = ut, as is 
readily checked. Condition (iv) and the network initial state set 5 = 0. Thus 

|* (t)> = cos^ |0) t |l) u |l> r |0) s + sin^ |l) t |1>„ |0) r |l) s (6) 

For if = | , one obtains ^ (^jj) = |l) t |l) u |0) r |l) s , namely the solution. 

The unitary transformation (6) brings the state of the network from satisfying only the 
input to satisfying both the input and the output constraints. It is obtained by "blindly" 
operating on divided parts of the network, but under A rs projection/s (the conquering factor) . 

The evolution is always unitary because an infinitesimal rotation of p s , under conditions 
(i) through (iv), yields a univocal (unitary) vector. We skip the lengthy but straightforward 
demonstration of this. 
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As a result of this process, A rs symmetries (or projectors) become constants of motion 
which commute with the network propagator at all times. They are also pairwise commuting, 
being applied to disjoint Hilbert spaces. However, the cause should not be confused with 
the effect. A rs projections shape or forge the unitary propagator with which they commute. 

If the network admits no solution, conditions (i) through (iv) make up an impossible 
system. Measuring the network final state — at t = j- — gives a non-solution. This is 
checkable in polynomial time and tells that the network is not satisfiable. 

If the network admits many solutions, the final state can be a linear combination thereof. 
Which one, depends on the network initial state through condition (iv). However measure- 
ment gives one solution (that it is a solution is checkable in polynomial time). 

It is clear from the above that Connections "cut" network complexity, inducing a divide- 
and-conquer strategy. This diakoptic approach would make NP-complete = P. However, 
we have applied reverse engineering until now. The A rs projections are just a nice-to-have 
feature. This raises the problem whether this feature is physical. 

V. INDUCED SYMMETRY 

A rs symmetry will be shown to be an epiphenomenon of fermionic antisymmetry in a 
special physical situation. This is generated by submitting a couple of identical fermions 
1 and 2 to a suitable Hamiltonian^ 12 !. We assume that each fermion has two compatible, 
binary degrees of freedom \ an d A. Just for the sake of visualization (things can remain 
more abstract), we can think that each fermion is a spin 1/2 particle which can occupy 
one of either two sites of a spatial lattice, x thus becomes the particle spin component a z 
(x — 0, 1 correspond to o z = down, up) and X = r,s the label of the site occupied by the 
particle. For example, |0) 1 11) 2 |r) 1 \s) 2 reads: a z of particle 1 down (0), a z of particle 2 up 
(1), site of particle 1 = r, site of particle 2 = s. 

The following is the list of the states which do not violate statistics; they make up the 
basis of the Hilbert space H\ x . States are represented in first and second quantization and, 
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when there is exactly one particle per site, in qubit notation (where o z j\ are the qubit 
eigenvalue/label) : 

|o> = 73 d°)i l 1 > 2 " l^i l°) 2 ) V)x \r) 2 = a l a\ r |0) , 



\b) = 75 d°>i l X ) 2 - l^i l°) 2 ) |s>i l«>2 = «£. a l |0> ; 

|c> = 75 l°)i |0>2 (lOl l«>2 " l*>i I0 2 ) = «Jr «L |0> = |0> r |0) s , 



|d) = ^ |l) x |1) 2 (Ir), | a > 2 - Is), |r> 2 ) = a\ r a\ s |0) = |l) r |l) s , 
|e) = |(|0) 1 |l) 2 + |l) 1 |0) 2 )(|r) 1 | S ) 2 -| S ) 1 |r) 2 ) = 
j= 2 (al a{ s + a{ r al) |0) = ^ (|0) r |1), + |l) r |0)J . 
l/) = H|0)il 1 ) 2 -|l)i|0) 2 )(|r) 1 | S ) 2 + | S ) 1 |r) 2 ) = 
j= 2 (al a{ s - a\ r aj.) |0> = ^ (|0> r |1), - |l) r |0)J . 

Creation/annihilation operators form the algebra: <a[,aH = {ai,aj} = 0, <a-,Oj} = Sij. 
Now we introduce the Hamiltonian H rs = E a \a) (a\ + E^ \b) (b\ + E c \c) (c\ + E^ \d) (d\ or, in 
second quantization 

H r s — ~{E a a-or ai r ao r a\ r + E\, a Qs a ls ao s ai s 
-\-E c a 0r a 0s ctor a os + E^ a lr a>i s a>ir a is) ■> 

with E a , Eb > E c , Ed > E discretely above 0. This leaves us with two degenerate ground 
eigenstates: 

\e) = ^ (|0) r |l) s + |l) r |0)J and |/> = ± (|0) r |l) s - |l) r |0>.) . 

Alternatively, their linear combinations |0) |1) and |1) |0) can be used as the two orthog- 
onal ground eigenstates. The generic ground state is thus: 

l^) = «|0) r |l) s + /3|l) r |0) s ,with |a| 2 + |/3| 2 = l, (7) 

which satisfies A rs symmetry. 

Let Ayi 1^) = \ (1 — Pyi) be the antisymmetrization projector. Due to the anticommu- 
tation relations: A 12 |0) r |1), = |0) r |l) s and A 12 |l) r |0) a = |l) r |0) s . Also, A 12 |0) r |0) s = 
|0) r |0) s and Ayi |l) r |l) s = |l) r |l) s , without forgetting that these are excited states. 
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The Connection can be implemented by suitably operating on the ground state (7). We 
assume that the initial ("symmetrical") state of the Connection is given by eq. (3): \^f (0)) = 
cosi?|0) r |l) a + sini?|l) r |0) s . Then transformation (2) p r (t) — Q r (u>t) p r (0) Q\ (ut) is ap- 
plied to qubit r, under continuous A rs projection. Let |\l/ (£)) be a free normalized vector 
of H\ x . The Connection state at time t is obtained by submitting \ty (t)) to the following 
mathematically simultaneous conditions, 
for all t: 

i) A 12 |tt(*)) = |tt(t)>; 

ii) Pr (t) = Tr s (|tf (t)) (* (t)|) = cos 2 (0 + <p) |0) r (0| r + sin 2 (t? + ^ |l> r (l| r ; 

iii) the distance between the vectors before and after reduction must be minimum, as 
specified in Section III; 

iv) the expected Connection energy: (£ (£)) = (\P (t)| iJ rs |\1/ (t)) , must be minimum. Since 
this minimum will always be zero, time ordering is irrelevant. 

It is readily seen that the solution of this system is still |\l/ (£)) of eq. (4): 
|*(t)) = cos(t? + ^)|0) r |l) a + sin(i? + ^)|l) r |0> a . 

Simultaneous satisfaction of (i), i.e. fermionic antisymmetry seen as projection, and (iv) 
(which is satisfied by (£ (t)) = 0) originates the projection A rs \ty (t)) — \^ (t)), as is readily 
seen. Therefore, if (£ (t)) = 0, namely if the operation on qubit r is performed adiabatically, 
we obtain the Connection. 

Since this computation is reversible 7 ' 8 ' , namely it does not dissipate free energy (the 
result of driving and shaping is a unitary evolution), in principle (£ (£)) can always be zero. 
This is of course an idealization, for the time being we are highlighting a possible, speculative 
way of dealing with NP-complete problems. 
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By the way we should note that the tensor products |0) r |0) s and |l) r |l) s that would be 
projected off since they violate particle statistics (Section II), are not the antisymmetrical 
excited states \c) and \d), which satisfy A&. They would be instead the symmetrical states 



1 
1 



|0)J0) S =-^|0> 1 |0> 2 (|r> 1 |s> 2 + |s> 1 |r> 2 ) 



|l) r |l) s =—\l} 1 \l) 2 (\r} 1 \s} 2 + \s) 1 \r} 2 ) 

The two kinds of states (antisymmetrical and symmetrical) have the same qubit notations 
and density matrices. We are of course in counterfactual reasoning; the important thing is 
that conditions (i) through (iv) yield the solution (4). 

Let us now address the problem of creating many Connections, namely an H rs Hamilto- 
nian per network wire r, s (fig. 2a). These H rs operate on disjoint pairs of qubits. Viewed 
as A rs projectors (which is the case when (£(£)) = 0), they are pairwise commuting. Still 
in the idealized case of adiabatic operation, the Connections operate independently of each 
other. 

VI. CONCLUSION 

The notion of applying a particle statistics symmetry (or projection) to divide the 
quantum whole into parts without clipping its richness — here computation speed- 
U p[9,io,n, among others] _ introduces an engineering (diakoptic) perspective in the design of 
quantum mechanisms. For the time being, the development of this idea remains at an ab- 
stract level. Finding model Hamiltonians which implement the Hermitean matrix of Section 
V could possibly be the next step. 

The interpretation of particle statistics symmetry as projection on a predetermined 
subspace is best modeled in a two-way (advanced and retarded in time) propagation 
scheme' 12 ' 13 ' 14 ]. 

This research has been supported by Elsag Bailey a Finmeccanica company. Thanks are 
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